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Abstract 

The correlated 3d sulphide BaVSa undergoes a sequence of three symmetry breaking transitions which are reflected 
in the temperature dependence of the magnetic susceptibility, and its anisotropy. We introduce a microscopic 
model based on the coexistence of wide band a\g and localized eg d-electrons, and give the classification of the 
order parameters under the double space group and time reversal symmetries. Allowing for the relativistic spin- 
orbit coupling, the d-shell multipoles acquire a mixed spin-orbital character. It follows that orbital ordering is 
accompanied by a change in the susceptibility anisotropy. 
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BaVSs can be regarded as a d-electron system with 
hidden order, or perhaps with several hidden orders. 
Its metal-insulator transition (MIT) has been known 
for over 25 years [1], and it was long supposed to be an 
example of a Mott transition without symmetry break- 
ing. No magnetic ordering accompanies the MIT and 
it was not discovered until 2002 that there is a struc- 
tural change: the unit cell is doubled along the trigonal 
axis [2]. The structural order parameter could be used 
to construct a Landau theory, and this would be suf- 
ficient to account for the fact that the MIT remains a 
second order transition in a wide range of pressure [3] . 
Nevertheless, we would like to consider BaVSs as a cor- 
related electron system, and find electronic (as well as 
structural) order parameters for its phases. A system- 
atic search for order parameters is based on symmetry 
classification. 

The high temperature symmetry group of BaVSa is 
Gs = PG^mmc^Tt where Tt stands for time reversal. 
The unit cell contains two V sites along the C3 axis, so 
the cell doubling mentioned previously means that the 
T < Tmit unit cell contains four V atoms [2,4]. 
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The symmetry is lowered from Gs by a sequence of 
three ordering transitions: a hexagonal-to-orthorombic 
transition at Ts = 250K, the MIT at Tmit = 69K 
(breaking translational symmetry) and finally, time re- 
versal invariance is broken at Tx ~ 30K. Certainly 
there is an antiferromagnetic ordering of the magnetic 
dipoles but it is probably accompanied by the ordering 
of composite spin-orbital octupoles. 

BaVSa is a 3d^ system. In the high-T hexagonal 
phase, the t2g level is split into a aig singlet, and an 
eg doublet of orbitals \a) and 1 6). It is a long-standing 
problem [1,6] whether only eg states are occupied, or 
both aig and eg levels are important. The current un- 
derstanding emerging from ARPES experiments [7] 
and DMFT calculations [8] is that aig and eg electrons 
are present in roughly equal proportions. The aig states 
form a wide band, while eg electrons carry localized 
(dipole, quadrupole, and arguably octupole) moments. 

The fact that the orbital degrees of freedom of the eg 
electrons play an important role is shown by the tem- 
perature dependence of the anisotropy of the magnetic 
susceptibility [6] . The source of anisotropy is that both 
t2g and trigonal eg electrons have unquenched orbital 
moment, and thus there is a first-order effect from the 
relativistic spin-orbit coupling (SOC). For vanadium 
compounds, the SOC is substantial (A « O.leV). Sus- 
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ceptibility anisotropy measured hy A — {xc — Xa)/Xa 
begins to deviate from at about 250K, reaches -.07 
at 70K, and +0.2 at low T (c is the trigonal axis, a one 
of the orthorombic axes) . The transitions at Tmit and 
Tx are marked by cusp-like anomalies of A. Published 
results allow a comparison between Xc and Xa only, but 
we predict further anisotropy within the a~b plane. 

The standard order parameter of the 250K transi- 
tion is purely structural: it is the change in the c/a ratio 
[5] . Lifting the degeneracy of the Cg level causes orbital 
polarization [11] which we measure by the expecta- 
tion value of the pseudospin Tz = (1/2) (X^s \o.S){aS\ — 
j6S)(6S'|), where S is the spin. The transverse com- 
ponents Tx, Ty of the orbital pseudospin are analo- 
gously defined. The time-reversal-odd component is 
Ty = (l/2)Lz, the orbital momentum along the trigo- 
nal axis. 

The evidence from the susceptibility anisotropy gives 
sufficient reason to assert that the orthorombic split- 
ting must be modest and therefore the orbital degree 
of freedom r remains important at all T's. 

We performed a symmetry analysis and preliminary 
mean field calculations for two cases: 1) considering a 
single V site as a unit, and 2) for the two-V-atom unit 
cell of the high-T structure of BaVSs. The details will 
be given elsewhere [10]; here we restrict ourselves to 
simple illustrative arguments. 

A single t2g electron in external fields can be de- 
scribed by 

iy = AS-F + A(l-L^)-gC),z-H-(L + 2S) (1) 

where the first term is the SOC, the second the trigo- 
nal field, the third is a quadrupolar mean field corre- 
sponding to the ordering of 

Oxz = —{LxLz + LzLx) (2) 

and the last is the Zeeman term. We arbitrarily chose 
one member of quadrupolar order parameter doublet 
{Oxz, Oyz}- We find that Oxz induces the mixed spin- 
orbital quadrupole 

0'xz^^{LxSz+LzSx) (3) 

The magnetic field effects are: 

(i) for H |[ (0,0, 1), in addition to Sz and Lz, also 
Sx and Lx are induced. 

(ii) for H II (1,0,0), in addition to Sx and Lx, also 
Sz and Lz are induced. 

(iii) H II (0, 1, 0) induces only Sy and Ly. 
Analogous results are found for other kinds of orbital 
order. In principle, we can hope to identify the type of 
orbital order from finding which non-diagonal elements 
of the susceptibility tensor have non-zero values. 

In contrast to the cubic eg, the trigonal Cg doublet 
carries both dipoles and octupoles. The argument given 



above is equivalent to stating that the spin-orbital oc- 
tupole TxSx — TzSy has the same symmetry as Sz ■ 
This may be relevant for the description of the nature 
of ferromagnetism of BaVSea [12]. 

Going over to the case of the two-atomic unit cell, 
we have to worry about the spatial distribution of Eg 
and a\g electrons. To a good approximation, there is 
one localized eg electron per two V atoms, i.e., one eg 
electron in each high-T unit cell containing the sites 
VI and V2. We developed a Kondo-necklace-type [9] 
spin-pseudospin model in which all ordering degrees of 
freedom are ascribed to eg electrons [10]. In addition 
to the spin (5*) and orbital pseudospin (r) degrees of 
freedom, we have a second pseudospin describing the 
position of the Eg electron in the unit cell {rjz = 1/2 
for VI, r]z = -1/2 for V2). The 8-dimensional local 
Hilbert space supports S'^ — 1 = 63 local order param- 
eters which we classified under the symmetry Gs = 
PQs.mmc'^Tt. 

We carried out mean field calculations in the S-t-t/ 
model, assuming a sequence of ferro-r and antiferro-r; 
ordering [10] . Alternating order of the rjz pseudospins 
amounts to the appearance of a four-V-atom supercell 
with the aig-eg-eg-aig character. Superexchange gives 
singlet binding at the eg-eg centers of the tetramers, 
explaining the observed spin gap [3] . Spin correlations 
are drastically reduced in the para-77 phase, though 
orbital order persists. 
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